Double symmetry breaking and 2D quantum phase diagram in spin-boson systems 
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The quantum ground state properties of two independent chains of spins (two-levels systems) 
interacting with the same bosonic field are theoretically investigated. Each chain is coupled to 
a different quadrature of the field, leading to two independent symmetry breakings for increasing 
values of the two spin-boson interaction constants Q,c and SI/. A phase diagram is provided in the 
plane (Ocfi/) with 4 different phases that can be characterized by the complex bosonic coherence of 
the ground states and can be manipulated via non-abelian Berry effects. In particular, when Q,c and 
0,1 are both larger than two critical values, the fundamental subspace has a four-fold degeneracy. 
Possible implementations in superconducting or atomic systems are discussed. 
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Spontaneous symmetry breaking is an important fea- 
ture of quantum phase transitions [l| and typically implies 
the appearance of degenerate ground states. An under- 
standing and control of the symmetries of a quantum sys- 
tem can lead to the realization of new quantum phases 
for fundamental studies and/or for quantum applications. 
One example is given by the Dicke model [2| describ- 
ing the coupling of a collection of independent tvifo-level 
systems coupled to a single-mode bosonic field. In the 
Dicke model, for a critical spin-boson coupling, a quan- 
tum phase transition occurs with the appearance of two- 
degenerate ground states Recently, such an Hamil- 
tonian has been experimentally implemented in an effec- 
tive way in a system where a pump-dressed Bose-Einstein 
condensate is embedded in an optical cavity [ij, iSj. In the 
case of the Dicke model, the Hamiltonian symmetry is the 
parity of the total number of excitation quanta in the sys- 
tem [j, 0] . The symmetry breaking across the transition 
means that the degenerate ground states (the so-called 
'superradiant' phase) are not eigenstates of the parity 
operator 11, contrary to the non-degenerate ground state 
(normal phase) obtained for sub-critical coupling. More 
generally, the spin-boson systems are currently attracting 
a great deal of attention, since they can be implemented 
with very controllable artificial systems such as dressed 
ultracold atoms in optical cavities and superconducting 
circuit QED0. 

In this Letter, we explore the properties of a spin- 
boson Hamiltonian describing two independents chains 
of two-levels systems interacting with the same bosonic 
field in a way related, but qualitatively different to the 
Dicke model. One chain is coupled to a quadrature of 
a bosonic field, while the other chain is coupled to the 
orthogonal quadrature of the same field via two indepen- 
dent spin-boson coupling constants, which we will call 
VLc and 57/ in the following. Like in the usual Dicke 
Hamiltonian, the parity operator H commutes with the 
system Hamiltonian and describes a symmetry of the sys- 
tem. However, we will show that H can be written as a 



product of 2 symmetry operators which can be indepen- 
dently and spontaneously broken if ilc" and 57/ are in- 
creased above the corresponding critical values VL'jj and 
■ We calculate the phase diagram and energy spec- 
trum of the elementary bosonic excitations of the system 
in the thermodynamical limit. A phase diagram with 4 
different zones in the (r^Ci ^i) plane is provided and char- 
acterized in terms of the macroscopic coherences. When 
both Q.C > ^^"^ ^ the ground state sub- 
space is four-fold degenerate, with complex bosonic co- 
herences. We show that adiabatic loops in the {flc^i) 
space provide non-abelian Berry geometric unitary trans- 
formations, which can be used to manipulate the four- 
fold degenerate ground states. Possible realizations in 
circuit QED and atomic systems are discussed. 

The 'double-chain' Dicke Hamiltonian under consider- 
ation is described by the following spin-boson Hamilto- 
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where is the bosonic creation operator with energy 
hujcav The total angular momentum operators and 
{JI and Jl) describe the ensemble of Nc {Nj) 
two- levels systems with transition frequency wji (wj). 
The first chain is coupled to the boson field quadrature 
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while the second is coupled to the quadra- 



ture i{a — flt), via the coupling constants flc and 51/ 
respectively. The collective atomic operators are de- 



fined for k e {/,C} as = Ef^tic^i' and 



and 



(l/2)(4 + J'l) = (l/2)Ez?=i(^'f + ^'^), with ai" 

^'^ ~^ Ti„..l: .„.-L_:.„, r„,. il,. l th pggudo-Spiu, SO 



a J: the usual Pauli matrices for the Ik 
that the angular commutation relations read [J^ , ] = 

±Sk,k'J^ ; [Jl, J^'] = 26k.k'J^ for k, k' e {/, C}. 

Let us start our analysis of such a 'double chain' 
Dicke model by describing the role of the parity oper- 
ator Hs in the standard 'single chain' Dicke Hamiltonian 
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Hs/h = iOcava'^a + io^J, + n/y/N{a+a^){J+ + J-) . The 
anti-resonant terms J_a and J+ prevent the conserva- 
tion of the total number of excitations + N/2 + a'^a 
(contrary to the Jaynes-Cumming or Tavis-Cummings 
models |8|), however the parity of the excitation number 
operator Ilg = exp(i7r[J2 -I- N/2 + a^a)) still commutes 
with the Hamiltonian[3|. In particular, XI^ transforms 
the operators in Hs as follows: 



n, 



Us : (a, J,) -^n,(a, J,)nt = (-a,~J,) 



(2) 



while Jz and a^a remain unchanged. Analogously, in the 
present 'double chain' Dicke model , the parity of the 
total number of excitations 11 is conserved. The corre- 
sponding symmetry operator is defined as: 

n = exp(i7r(a^a + jf + Nc/2 + jI + Ni/2)). (3) 

It is apparent that [H, 11] = because: 

n : (a, Jf , 4) n(a, jf , Ji)Il^ = (-a, - jf , - 4). (4) 

The novelty given by the coupling to the two differ- 
ent quadratures is the possibility to express the parity 
operator as a product 11 = 7/ o Tc where 7} and Tc 
are independently conserved and are defined through the 
following transformations: 

(a + a),i{a - a)), jf , J^) ^ (a + a\-i{a - a)), , ~ J^), 
(a + a\ i{a - a^), jf , 4) 5 {-a - a\i{a ~ a^), - jf , 4) 

where a^a, Jl and remain unchanged by those two 
transformations. 7} and 7c are one the dual of each 
other: each one flips the sign of the corresponding 
quadrature of the field, while keeping unchanged the 
other one. There is a simple geometric interpretation 
of those symmetries in the complex quadrature plane 
(i?e((a)), Im{{a))). H is represented by a 7r-rotation with 
respect to the origin in such plane, Tc is a mirror reflec- 
tion with respect to the imaginary axis, while 7} is a 
mirror reflection with respect to the real axis. As we 
will show, the symmetries associated to Tc and Ti break 
down when the corresponding coupling constants Hp and 
ri/, are increased across some quantum critical lines. 

To explore the properties of such an Hamiltonian, it is 
convenient to use the Holstcin-Primakoff transformation 
which allows us to represent the angular momentum 
operators in terms of 2 independent bosonic modes: = 

biiN, bib.r/^ , ^ {N, 
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FIG. 1; Phase diagram of the 'double chain' Dicke Hamilto- 
nian in Eq. exhibiting four different phases in the two- 
dimensional parameter {Q,c,0.i) space, where fie {^i) is the 
spin-boson coupling for the first (second) chain of pseudospins 
(two- level systems). For each of the 4 phases, a sketch of 
the complex boson coherences (a) = ±.^70 + ^i^/li for the 
ground states is depicted in the complex plane (blue points). 
In particular, for Q.c > i^c ^^'^ > the phase is 

characterized by 4 degenerate vacua, with 4 different com- 
plex bosonic field coherences. Such coherences are mutually 
related via the symmetry operations Tc and Ti defined in the 
text. 



Note that the energy of the bosonic excitations in the 
normal phase can be obtained in the thermodynamical 
limit (Nk — >■ 00) by neglecting the terms proportional 
to {1/Nk) (for k = C,I) in the last Hamiltonian which 
therefore becomes quadratic in the three bosonic modes 
a, hi and be- The frequency of the lower energy mode 
can vanish for flc = ^c ~ {^/'^)\/^cav^c ^^'^ ~ 



(1/2) ^yu}cav(^^I, hence we have here two straight 
critical lines (one horizontal and one vertical) depicted 
in Fig. m 

To determine the phase diagram, we express the op- 



erators as bi 



blbkY'^k and 4- 



blpk — Nk/2 for k £ {/, C}. In this representation, the 
Hamiltonian reads: 



d\ - VPk {k = C, /) and c) + 

/"1c + where 7c 7/, /3c and /?/ are real and such 
that "fc ^ Pc ^ Nc and 7/ ~ /?/ ^ TV/ Q. By making 
the corresponding substitution in Eql5] and developing 
the square roots (1 - (4 - VPk){dk - ^/^)/Nkf/^, we 
get an Hamiltonian with constant, linear and quadratic 
functions of the bosonic operators c, dc , dj and their her- 
mitian conjugates. By setting the linear terms to zero Q 
we flnd the solutions for the spontaneous mean fields: 



_ ncVNc{i~^l) , . ni^Ni{i~fi-j) 

7c + iVT/ = ec — + 



H/h= ujcava^a + a;p(&|^&c 



Nc. 



u"i{b\b 



Ni, 



+ncia + a^){bl^l -bl.be/Nc + ^l-b^bc/Nc be} 



+ini{a - a^){blJl - b\bi/Ni + V 1 - b]bi/Ni bi}. (5)tion fik 



= ec^^(l - Ac) ; = e/^^(l - M/). (6) 
To write the results in a compact way, note that by dcfini- 

n 



lif rifc < f^^''; Afc 



4n? 



if Ofc > , while 
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FIG. 2: Frequencies of the 3 branches of bosonic excita- 
tions (called lower, middle and upper) lji < u),n < Cju for the 
resonant case cj? = uj^ = ojcav and in units of ujcav 3D sur- 
face plot: lower eigenmode ui as a function of Qc and fi/. 
The other plots depict the three excitation frequencies versus 
^/i^Q + Q,'j/ujcav for particular values of the polar angle in the 
parameter plane, namely 6 = 2n /5, d = n /4 (up) and 9 = Tr/7 
(down). For 9 ^ n/A [7r/2], two different critical points mean 
two consecutive symmetry breakings. For 6 = 7r/4 [tt], the 
two symmetry breakings occur at the same time, implying a 
single quantum critical point. 



Efc = ±1. These equations predict the existence of 4 dif- 
ferent phases. For ilc < and ftj < iiy , the normal 
phase is characterized by a non-degenerate vacuum with 
no macroscopic coherence {-jc = 7/ = /?c = /?/ = 0). 
For flc > ^^'^ < ^ there is a superradiant 
phase (like in the Dicke model) with a twice degenerate 
vacuum with a real boson coherence (a) and a macro- 
scopic spin polarization (along the .x-direction) for the 
chain labeled by C, which is coupled to the quadrature 
(a -(- flt). The third phase corresponds to ilc < and 
51/ > riY- it is also superradiant with a twice degenerate 
vacuum with imaginary boson coherence and a macro- 
scopic polarization of the chain of pseudo-spins labeled 
by /, which is coupled to the quadrature i{a — a^). Fi- 
nally, when Qc > ^^e phase is dou- 
bly superradiant with a fundamental subspace that is 4 
times degenerate, with complex bosonic coherence and a 
macroscopic polarization for both the two independent 
spin chains. Such a phase diagram is summarized in Fig. 

m 

In order to determine the energies of the bosonic exci- 
tations for the different phases, i.e. for any value of 51/ 
and He, one has to diagonalize a bosonic Hamiltonian 
quadratic in the fields c, dc, di and their hermitian con- 
jugates. There are 3 branches of excitations whose posi- 
tive frequencies (uji < LOm < Cby) are obtained by diago- 

V Q 



nalizing a 6 x 6 Bogoliubov matrix Jv{ 



where the 3x3 hermitian submatrix V and the 3x3 
symmetric submatrix Q read: 



V = 




51f 



iil/ 





+ 2bc 

















— iil/ 


Oc 


-2Dc 
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-2Di 



2D, 



(7) 



(8) 



with 



^"(1 + Afc)/(2/ifc), flk = \/2flfc/ife/Vl + Mfe 
and = {cjO(3 + /2fc)(l-/ife)}/(8Afc + 8/i2) {k = C,I). 
The values of the excitation frequencies as a function of 
(lie, 51/) are shown in Fig. [2] The behavior of a); at the 
quantum critical points is analogous to the one in the 
standard single-chain Dicke model 0]. The double criti- 
cal point {Qq, flf ) behaves differently as the frequency 
vanishes linearly with the coupling. The general expres- 
sion (valid for any phase) for the ground state energy 
reads: 



Er 



^ ^{iV,(l+M^) + (l-A.)}. 

k=C,I 



(9) 



It is also interesting to show a very simple asymptotic 
expression of the 4 degenerate vacua of the 'double su- 
perradiant phase', in the finite size case, which we have 
derived in the limit Nk ^ (flfc/wcau)^ {k = C, /): 




which is a product of a coherent state for the boson 
part times two ferromagnetic states polarized in the 
pseudospin cc-direction for the two chains, namely | ± 
Nc/2)^'^\±Ni/2)^. Those states satisfy I ±iVfc/2)^ = 
±(A^fc/2)| ± Nk/2)x {k — C,I), meaning that each atom 
is maximally polarized. For a finite number of spins, 
the exact degeneracy is lifted, the first 4 eigenstates be- 
ing linear superpositions of the 4 states \G±±). The 
3 energy splittings which appear exponentially decrease 
with the size and/or the coupling constants, being either 
~ exp(-2Afc^^?7/w?aJ - cxp(-_2iV/5l2/u;2^J in anal- 
ogy to the standard Dicke model ll| , implying that for 
many applications they can be neglected. 

An interesting property of the present system is the 
possibility of manipulating the four-fold ground space via 
geometric Berry effects. If we vary 51c and 51/ to form a 
loop trajectory C in the corresponding two-dimensional 
parameter space (as shown in Fig. [3]) in an adiabatic 
fashion (slowly with respect to 1/ujcav), then one creates 
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FIG. 3: Left: Closed loop C enclosing an area A in the 2D pa- 
rameters plane (fie, If the trajectory occurs in the part 
of the phase diagram with the four-fold degenerate ground 
state, then the corresponding adiabatic transformation pro- 
duces a non-abelian Berry unitary transformation in the basis 
{\G++),\G+-),\G-+),\G—)}. ' 



an non-abelian Berry unitary transformation [I3 U{A) 
in the basis {\G++), |G+_), |G_+), \G )}. Namely: 



U{A) = Vt 



(11) 



where the geometric ang le is ~ 2A^/NcNi/ujI^^, with 
A the loop area in the parameter plane {^c^ ^i) (see Fig. 
E]) and (g) = dmg(l, -1, -1, 1). The path ordering 
operator V had been be taken away because the Wilczek- 
Zee connections [1^ {G ±±\-i^\G ±±) arc almost diagonal 
since exp(-2Affc^2|/w^^„) < 1 (for k = /,C). Note that 
if one sees the two-chain system as a two-qubit system, 
then such non-abelian transformation is equivalent to a 
conditional two-qubits quantum gate. We would like to 
point out that we have successfully tested our analytical 
results with finite-size exact diagonalizations. 

The double chain Dicke model may be realized in cir- 
cuit QED systems, by considering, for example, a chain 
of Josephson two-level atoms capacitivcly coupled to a 
transmission line resonator and considering another chain 
of Josephson atoms inductively coupled to the same res- 
onator. The 'capacitive' and 'inductive' schemes are 
two different ways of coupling a two-level system in a 
transmission line resonator in circuit QED[13|. The ca- 
pacitive coupling connects the charges of a Josephson 
two-level system to the resonator voltage quantum field 
Vr oc {a + a^) 0. The inductive coupling, which is dual 
of the latter , provides an interaction between the flux 
of the Josephson atom and the current of the resonator 
which is proportional to the other quadrature of the field, 
i.e. Ir oc i{a — a^), because the resonator current is pro- 
portional to the time derivative of the resonator voltage. 
In particular, a chain of Ni identical fluxonium Joseph- 
son atoms [IJ], inductively coupled to the first bosonic 
mode of the resonator can undergo a spin-boson inter- 
action of the form Hi/h = 2il7////V7(a - a^)Ji [13. 
The coupling Hamiltonian of a collection of Nq identical 
Cooper pair boxes [l5|. interacting with the first bosonic 
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FIG. 4: A possible realization of the Hamiltonian [T] in cir- 
cuit QED. One chain of Cooper pair boxes (blue) capacitively 
coupled to the voltage bosonic field of a transmission line res- 
onator and one chain of fluxoniums (green) inductively cou- 
pled to the flux bosonic field of the resonator. Because of the 
boundary conditions 0], the first mode profile of the voltage 
(resp. current) field has its antinodes at the edges (resp. cen- 
ter) of the resonator, as shown above in solid (resp. dashed) 
line. The capacitive coupling provides a coupling to a quadra- 
ture of the resonator bosonic field, while the inductive cou- 
pling provides the coupling to the other quadrature. 



mode of the resonator, with a gate capacitance Cg and a 
static gate voltage Vg such that Ug = CaVa/{2e) = 1/2 



can be written as Hc/h = 2f]c/V^(a + at) jf jl6l[l7|. 
One way to combine these two kinds of couplings and 
get the Hamiltonian [1] is sketched in Fig. SI Note that in 
such a context, the transformation 7/ might be defined 
as the time reversal symmetry since it lets unchanged the 
resonator voltage and the qubit charge operator while it 
flips the resonator current and the qubit flux. And Tc 
will be then the product of the parity transformation IT 
and the time reversal symmetry 7/. Finally, this kind 
of Hamiltonians could also be implemented with atomic 
Bose-Einstein condensates dressed by judiciously tailored 
laser pump beams in an optical cavity [l8| in a configu- 
ration related to the recent observation of the dynamic 
Dicke quantum phase transition^, 



In conclusion, a spin-boson system with a quantum 
phase transition characterized by two independent sym- 
metry breakings has been presented. The full phase dia- 
gram of such 'double chain' Dicke model has been solved, 
showing the possibility to obtain four-fold vacuum de- 
generacy associated to a 'doubly superradiant phase'. 
We have shown how to manipulate the ground states 
via non-abelian Berry methods. Our work provides a 
new paradigm in the very active field of spin-boson sys- 
tems, with possible applications in circuit QED and cav- 
ity QED with superfiuid atoms. We wish to thank K. Le 
Hur, J. P. Gazeau, H. Ritsch, D. Hagenmiiller for fruitful 
discussions. C. C. is member of Institut Universitaire de 
France. 
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